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OO 

We investigate the power-law tails in the evolution of a charged massless scalar field around a 
fixed background of dilaton black hole. Using both analytical and numerical methods we find the 
inverse power-law relaxation of charged fields at future timelike infinity, future null infinity and along 
the outer horizon of the considered black hole. We invisaged that a charged hair decayed slower 
than neutral ones. The oscillatory inverse power-law along the outer horizon of dilaton black hole 
is of a great importance for a mass inflation scenario along the Cauchy horizon of a dynamically 
formed dilaton black hole. 
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I. INTRODUCTION 



The late-time evolution of various fields outside a collapsing star plays an important role in major aspects of black 
1 hole physics, i.e., in the Wheeler's no-hair theorem Q and in mass inflation scenario In Ref. || the neutral 
external perturbations were studied and it was found that the late-time behaviour for a fixed r is dominated by the 
factor t~( 2l+3 \ for each multipole moment I. Next, Gundlach et al [Q studied the behaviour of neutral perturbations 
along null infinity and along future event horizon. Their conclusion was that along null infinity the decay is according 
the power law of the form u~^ l+2 \ where u is the outgoing Eddington-Filkcnstein (ED) null coordinate. On the other 
hand, the neutral perturbations along the event horizon behave like w _ ^ +3 \ where v is the ingoing ED coordinate. 
Bicak || analysed the scalar field perturbations on Reissner-Nordstrom (RN) background and found for | Q |< M 
the relation t~^ 2l+2 \ while for | Q |= M the late-time behaviour for a fixed r was governed by the rule t~( l+2 \ 

The late-time behaviours of a charged massless scalar field in RN spacetime were studied analytically and confirmed 
by the direct numerical studies of the fully nonlinear gravitational collapse in |^ ^| . Among all, the conclusion was 
that a charged hair decayed slower than a neutral one, i.e., the charged scalar hair outside a charged black hole was 
dominated by a i~( 2i + 2 ) behaviour. 

The problem of the late-time tails in gravitational collapse of a self interacting massive field were studied in Q (see 
also |fi0|-p2[). At late-times the decays of a self-interacting hair is slower than any power law. 
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In our work we shall discuss the asymptotic evolution of a massless charged scalar field in the background of dilaton 
black hole. Focusing on no-hair theorem we shall be interested in the dynamical mechanism by which the charged 
hair is radiated away. In Sec. II we start with the outline of the system under consideration and we analitycally study 
the late-time evolution of charged scalar perturbations along timelike infinity, future null infinity and the black hole 
outer horizon. Then, in Sec. Ill we treated the problem numerically. Sec. IV outlines our conclusions and remarks. 



II. THE EINSTEIN-MAXWELL-DILATON EQUATIONS 



In this section we shall analytically study the evolution of a massless charged scalar field tp around a fixed background 
of electrically charged dilaton black hole. The wave Eq. for a massless charged scalar field is given by [Ol 



V^V^V - ieA a g al3 (2Vpip - ieAptp) - ieV^A^ = 0, 



(1) 



where e is a constant value and Ap is the gauge potential. 

The metric of the external gravitational field will be given by the static, spherically symmetric solution of Eqs. of 
motion derived from the low-energy string action (see e.g. 111]). The action has the form as follows: 



S = J d A x [R - 2(V0) 2 - e^F 2 ] 



(2) 



where <j> is the dilatonic field and F a p — 2V\ a Am. The metric of the electrically charged dilaton black hole is as 
follows: 

dr 2 



1-^ 



df z 



r I 



r -£)<*" 



sin 



(3) 



The event horizon is located at r + = 2M. For the case of r_ = ^ we have another singularity but it can be ignored 
because of the fact that r_ < r + . The dilaton field is given by e 2 ^ = e~ 2 ^° (l — —J, where 4>q is the dilaton's value 



at r — > oo. The mass M and the charge Q are related by the relation Q" = — - — < 



2 _ r + r- 2 <t>n 
2 



Defining the tortoise coordinate y, as dy 



2M 



dr 



one can rewrite the metric (0) in the form 



ds = ( 1 - — ) H* + d v J + M r ~ ) ( dff 



Q 2 



'I- 



(4) 



For the spherical background each of the multipole of perturbation field evolves separetly so for the scalar field in the 
form ip = J2i m Vm(t' r )^ m (^ 4 , )/R( r )i one nas the following Eqs. of motion for each multipole moment 



ry jtt - 2ieA t rj. t - T), yy + Vr\ = 0, 



where 



V 



1(1 + 1) R" 



R 2 



R 



2M 
r 



2MR' 
r 2 R 



1-^1 



e 2 A 2 . 



r I 



(5) 



(6) 



By R we denoted R 



M 



and ' is the derivative with respect to the r-coordinate. As in Ref. Q in order 

Q 



to get rid of the physically unimportant quantity which appear in A t = $ — ^f, one can define the auxilary field 



^ — ie<£t 



i]. Then Eq.(|^) may be written as 
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Q - 

ip t u + 2ze— ip tt - ip :Vy + Vip = 0, 



where 



V = 



1(1 + 1) R^_ 
R 2 + ~R 



1 - 



2 ") 



2MR! 
r 2 R 



1 - 



2M 



e 2 Q 2 



One can assume that the general solution of Eq. (0) is of the form 

l 

# = ^A k R- k \ e - ieQlnr G(u) {l - k) + (-i) k e ^Q lnr F(v) {l ' k) 



k=Q 



(7) 



(8) 



(9) 



OO 

[B k (R)G{uf- k -V + C k F(uf- k -V 



k=0 

where F and G are arbitrary functions of an retarded time coordinate u = t — y and an advanced time coordinate 
v = t + y. The superindices on F and G have different meaning depending on negativity or positivity of them. The 
positive superindices indicate the number of times the function is differentiated, while the negative ones are to be 
interpreted as integrals ||. The first sum in ijj represents the primary waves in the wavefront while the second sum 
depicts the backscattered waves. The first sum in the above expression is the zeroth order solution. 
The recursion relation for B k (R) is as follows: 



2X 2 B k (R)' 



^B k (R)-\ 2 [X 2 B k ^(Ry}' + X 
V [A k (R)R- k r- ieQ 



A k (R)(kR' + ieQ) 



1 



(10) 



\ieQ(R-r\ 2 ) - \ 2 rR'(k + l) 
B k _ 1 (R)]+2A k+1 (R)^-^ „„/, V ± = 0. 



Jj>fc+l r ieQ+l 



where A = (1 — We expand the B k (R) coefficient in the form 

dk b k 



B k (R) 



J^k+lj,ieQ J^k+lj.ieQ+1 

where as in (|] ak = a k (l, eQ), b k = b k (l, eQ). 

Using Eqs.(^0|) and ( pi] ) the lowest order coefficients yield 



(11) 



a k = -ieQA k (R) 



2k + 1 



[l + 0(eQ)} 



(12) 



2(k + 1) 

As was pointed out after starting of collapsing of the star surface it approached very fast the speed of light. Thus its 
word line is asymptotic to an ingoing null geodesies v — vq. The time dilatation between static frames and infalling 
frames caused that the variation of the field ip on the star's surface is asymptotically redshifted. This means that 
ip tU will be small at late times. Then it is well justified to make an assumption that after some retarded time u = u\ 
variations of the charged field ip on v — Vq can be neglected. 

In the first stage of evolution we shall consider the scattering in the region u < u < u i- One has taken the variation 
of the field on v — vo to be neglected after u\. Thus, we have no outgoing radiation for u > U\, except from 
backscattering. This assumption about no late outgoing waves is equivalent to taking G(u\) = 0. Summing it all up, 
for large r at u = u\ the dominant term in ip (the dominant backscatter of the primary waves) is as follows: 



4>(u = u 1 ,r) = ai-^G { 



(13) 
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In the next stage of the consideration we take a closer look at the asymptotic evolution of the charged scalar field. 
We shall take into consideration the region where y ^> M, Q and for which u > u\ . In this region the leading order 
effect on the propagation of ip, has the centrifugal term in the potential V, namely the evolution is dominated by a 
neutral flat spacetime terms. Thus one has 

ttt ~ VVr + = 0. (14) 

As in RN case in this region one can replace r by y, because of the fact that ip,yy — ^,rr + 0(M). The solution to 
Eq.(|l4|) can be written in the form 

k 



i> = J2A k y- k g«- k Hu) + (-l) k f {l - k Hv) . (15) 



Matching ( pL 5[ ) with initial data on u = u\ Eq.(|13[), one has f(v) = Fqv~ , where Fq is of the same form as in RN 
case The standard procedure of expansion jjj, for late times t 3> y for g(u) = X^^Lo ^~~r ff^Wz/" an d similary 
f( v ) — 12^=0 enables one to rewrite ( |l5| ) as follows: 

oo 

$ = J2 K ^ n [/ ( ' +n) (i) + (-l) n 9 {l+n) (t)] , (16) 

n=-l 

where K n are given in Ref. || as in the neutral case. 

Using the boundary conditions for small r (see || and Q for the details), one finds that the late time behaviour of 
the scalar field at timelike infinity i + is 

4> ~ 2K l+1 y l+1 f < - 2l+1 \t) = -2K l+1 F (2l + l)\r 2 ^y l+1 + O(eQ). (17) 

On a future null infinity scri + one has 

ip{v > u,u) ~ A g( l) (u) ~ -F Z!ir (l+1) . (18) 

At the black hole outer horizon H + (as y — > — oo) the Eq. of motion can be approximated by 

^, t + ^,-^yy-^ = Q, (19) 

with the general solution of the form 

ip = e~™" [a(u) + j(v)\ = e~~ *y(v). (20) 

We put a(u) = because of the condition of neglecting the variation of the scalar field on v = v Q . To join the solution 
for y <c — M with y ^> M for i + , one uses the ansatz |^,^| ip ~ "ipstatt^ 2 ^ +1 ^ for the solution in region y <C — M 
and i 3>| y |, i.e., we assume that for the whole range of y the solution has the same late-time dependence as for 
y M. This enables us to match the solution (|2(]) with that for i+. One concludes that the late-time behaviour on 
the horizon H + is 

$(« -> oo, v) = r e^w~ 2( ' +1) , (21) 
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where the constant Tq is detrnhned by the condition that there is a static solution fulfilling /im J( ^ 00 V'/i 



Tq and 



lirriy-^oai) finatat = ~ 2Kf l+1 Fq(1 + 2l)\r l+1 . As was remarked in ||f|], the coefficient To and all other aspects of the 
late-time behaviour are governed by the initial backscattering. 

As one can see the behaviour of the charged scalar field on a fixed dilatonic background at timelike null infinity i + 
and the future null infinity scri + is the same as in the RN black hole case. However, on the external horizon H + we 
have some changes, but the general form of the time dependence is maintained. 

III. NUMERICAL METHOD AND RESULTS 

In this section we shall numerically study Eqs.(0) and (||). First we transform Eq.(^) into (u,v) coordinates. This 
yields: 



4$ )W) + 2ie—(ii, u + ip, v ) + VtJj = 0. 



(22) 



Then we divide this Eq. into a set of coupled equations for real, ijijn, and imaginary, ipj, parts of the field. 

4ii K , uv - 2e-(Vv + + Vifa. = 0. (23) 
r 

4^ <uv + 2e-(^ R , u + i> R , v ) + V$i = 0. (24) 
r 

To solve these equations numerically we establish uniformally spaced grid («,,«,-), i — 0, ...,n, with Ui = Vi = ih and 
Au = Av = h. Next, we construct implicit difference scheme: 



(25) 
(26) 



where if ip s denotes the value of the considered field at point with coordinates (uo,vq), ip w , w, w are the values of 
the field at (uq + h, vq), (uq, vq + h), (uq + h, vq + h), respectively (see Fig. [I]). 
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FIG. 1. Schemat of computational grid. Initial conditions are specified on u = uq axis, and boundary conditions on v = vo 
axis. Given the values of the field at W,S and E points we can calculate the field at N. Starting from lower corner of the grid 
calculation proceeds to upper corner row by row. 



Thus given the initial values of the field on the u = uq axis and imposing boundary conditions on the v = v o axis 
one may integrate Eqs. (^5|) and (^) row by row. As an initial condition we choose the Gaussian impulse of the form 

{v-v*) 2 ~ 



ip = (1 + i) exp 



(27) 



In all calculations presented here we use v* = 100 and a = 400. The size of the grid ranges from to to 10 4 or 
4 x 10 4 in both u and v and usually h equals at least 0.1. Because of the scale invariance of the problem we choose 
2M = 1. Unless otherwise noted we use Q — 0.45 and e = 0.01. 

In Fig. ^ we present time evolution of the field on three different hypersurfaces: future timelike infinity i + , future 
null infinity scri + , and the future horizon H + . In the numerical experiments these lines are approximated by the 
values of the field ip(y = 400, t), ip(v — 4000, m), and %j)(u — 2 x 10 4 ,t>) respectively. Calculations are done for I = 0. 
After a short period of a quasinormal ringing the dominance of power-law tails becomes prominent. The solid lines 
indicate power law exponents —1 and —2 expected from the theoretical predictions of Sec. O. As can be seen the 
agreement between these predictions and the experiment is excellent. The highest discrepancy for scrit + results from 
the fact that theoretical assumption w > u is not well maintained at the end of the line. 




FIG. 2. Time evolution of the scalar field from the Gaussian initial data. i+ approximated by tp(y = 400, t) represents future 
timelike inifinity while scrit+ (4>(v = 4 x 10 4 ,w)) and H+ (4>(u = 2 x 10 4 ,u)) represent future null infinity and the future 
horizon, respectively. Numerically calculated slopes for these curves are —2.04 for i+, —1.13 for scrit+, and —1.99 for H+. 
Thin, solid lines have slopes equal to power-law exponents expected from theoretical predictions. 

In Fig. |^ we examine the dependence of power-law exponents on the Z-number. We studied the behaviour of the field 
on future timelike infinity i+ for I = 0, 1, 2. As in the previous example the thin solid lines indicate the theoretically 
predicted values of power-law exponents from Eq. ([l7|). Once again the late time the power-law tails of the field 
converge to those theoretically predicted. 



G 



10 3 10* 
time, t 

FIG. 3. Field evolution on future timelike infinity, i+, for different values of the parameter I — 0, 1, 2. Curves have slopes of 
— 2.04, —4.07, and —6.06, respectively. Thin solid lines indicate power- law exponents predicted from analytical investigations. 

In order to check the late time behaviour of the field depends on the charge Q we performed the calculation for two 
values of (Q/M) 2 ratios 0.1 and 0.9. For the Gaussian initial data with I = 0. The results are presented in Fig. |J. 
Once again the theoretical predictions arc in agreement with the numerical simulation. Fig. |] envisages the evidence 
for the existence of power-law tails on the outer horizon of dilaton black hole for generic perturbations being of a great 
importance for the mass inflation scenario. For the first time the existence of power-law tails on the outer horizon of 
RN black hole was revealed in Q studing a scalar field perturbations and also confirmed for the case of a massless 
charged scalar fields on RN background |^,|| . 




FIG. 4. Field evolution on future timelike infinity for different values of (Q/M) 2 . Power-law exponents of the late-time tails 
are -2.07 for (Q/M) 2 = 0.1 and -1.98 for (Q/M) 2 = 0.9. 
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IV. CONCLUSIONS 



In our work we studied the asymptotic late-time behaviour of a massless charged scalar held in the background of 
dilaton black hole being the spherical solution of the so-called low-energy string theory. We confined the existence of 
oscillatory inverse power-law tails in a collapsing spacetime along asymptotic regions of future timelike infinity, future 
null infinity and along the outer horizon of the considered black hole. 

On a null grid (u, v) we integrated numerically the linearized charged scalar field Eqs. and confirmed the existence 
of the power-law tails. One also treated the dependence on the Z-number, they converged to the results theoretically 
predicted. The numerical results confirmed our analytical conjecture that a charged hair decayed slower than a neutral 
one. The same statement was also revealed in RN case §,§. The oscillatory inverse power-law tails along the outer 
horizon of dilaton black hole suggests the occurence of mass inflation along the Cauchy horizon of a dynamically 
formed dilaton black hole. It will be interesting to study the fully nonlinear gravitational collapse using a different 
scheme based on double null coordinates [jl5|] which allows one to begin with regular initial conditions and continue 
the evolution all the way inside black hole. We hope to return to this problem elsewhere. Another problem for 
investigations is the problem of late-time behaviour of self-interacting scalar fields in the spacetime of dilaton black 
hole and the higher order string corrections modifying the initial action. 
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